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Abstract 

We point out that, in some models of small-x hard processes, the transition 
to the Pomeron regime occurs through a sudden tunneling effect, rather 
than a slow diffusion process. We explain the basis for such a feature and we 
illustrate it for the BFKL equation with running coupling by gluon rapidity 
versus scale correlation plots. 
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1 Introduction 



In most situations in QCD, the question of whether a given observable can be calculated 
perturbatively reduces to two issues: firstly whether all the scales in the problem are 
large compared to Aqcd and secondly whether the observable is infrared and collinear 
safe. 

One notable exception is the limit of high-energy scattering of two hard probes, for 
example quarkonia, or virtual photons. Fixed order analysis suggests that the scale of 
a s is simply determined by the inverse transverse size of the objects being scattered, 
independent of the centre of mass energy. But at very high centre-of-mass energies, 
y/s, it becomes necessary to consider diagrams at all orders in the coupling. This 
resummation was initially considered for fixed coupling |l|], but subsequently interest 
developed in its extension to the (more physical) running coupling case 0, |3|. One 
of the conclusions of these studies was that regardless of the hardness of the probes, 
there always exists a centre of mass energy beyond which the problem becomes entirely 
non-perturbative, see for example @, || . 

Other next-to-leading corrections 0, which however must be resummed, as is done 
for example in f7|, §, || or in other related fH| |ll|] or less related [|l^, [13| approaches, 



do not qualitatively change this picture. The transition to the non-perturbative regime 
takes place also in that case, although the diffusion effects can be different, due to the 
lower value of the intercept and the diffusion coefficient. This problem is now under 



detailed investigation [14 



For a long time, the view of the transition to the non-perturbative regime was based 
on the idea of diffusion. At lowest order in a s the dominant diagram at high energies is 
the exchange of a gluon between the two probes. The scale for a s is determined by the 
transverse momentum of the exchanged gluon, k 2 , typically of order of transverse scale 
of the probes, Q 2 (in what follows we will generally refer to t = lnQ 2 /A 2 ). At higher 
orders one needs to consider emissions of final-state gluons from the exchanged gluon 
— these emissions modify the transverse scale of the exchanged gluon, which undergoes 
a random walk in t' = Ink 2 / A 2 . In [15, 16] this diffusion was illustrated graphically by 
plotting the median trajectory in t' as a function of rapidity (the distance along the 
chain) as well as the dispersion around that median trajectory. A schematic version of 
such a plot is shown in figure [l|, with the dark (blue) band representing the extent of 
the diffusion. Its characteristic shape led to it being named a cigar fT6|j . 

While at its ends the cigar is restricted to being narrow (by the presence of the 
probes), in the middle its width grows as a function of the centre of mass energy: 
At' ~ \j6t s Y with Y = In s/Q 2 . For the problem to remain perturbative, typical 
trajectories must remain well outside the non-perturbative infrared region and this 
leads to the requirement that the width of the cigar be less than t, i.e. y/a s Y <C t. 
Using a s (t) = l/bt, with b = 11/12 the coefficient of the 1-loop ft function, this reduces 
to the condition Y <C bt 3 . This limit was originally derived (with some additional 
numerical factors) in |T7|, [T3], ^tj. 



This picture of diffusion has its roots in the fixed-coupling approximation. It is 
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Figure 1: Above: a chain of emissions stretching in rapidity (0 < y < Y = In s/Q 2 ) 
between the two hard probes at scale t = \nQ 2 /A 2 . Below a schematic representation 
of typical transverse scales at different stages along the chain: the white line illustrates 
a typical random walk in t' (Ink 2 ) as one goes along the chain, while the dark (blue) 
'cigar' represents the average limits of this diffusion. 
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tempting to consider the full running- coupling case as a small variation of the fixed- 
coupling situation. One then expects the differences to be that the centre of the cigar 
dips to t' values below t, and that the cigar profile in t' develops some asymmetry with 
respect to its centre (the cigar becomes a 'banana', figure 0a). This approach to the 
running-coupling case as a perturbation of the fixed-coupling case leads also to well- 
defined analytical predictions [18|, [19], ^(| [Hj]. In particular, the usual fixed-coupling 
exponential growth of the cross section, log a ~ a s Y, is supplemented by an extra term 
proportional to b 2 a^Y 3 . This can be interpreted as due to the centre of the cigar being 
at if values of t — O (ba 3 Y 2 ), which for Y of order bt 2 implies that the cigar reaches the 
non-perturbative domain. Indeed calculations of the 'purely perturbative' component 
of the cross section show unphysical (oscillating) behaviour for Y > bt 2 , which can be 
interpreted as meaning that the perturbative component is not a physically well-defined 

22l. 



quantity in that region [21 




Figure 2: Analogues of fig. [I] which show the possible effects of the running coupling on 
the shape of the 'cigar'. 



Recently however, it was discovered in the context of a simplified, collinear model for 
BFKL dynamics [^3 , that the running of the coupling can introduce further qualitative 
differences. Rather than the width of the cigar gradually increasing with Y until it 
reaches the non-perturbative region at Y ~ bt 3 , or the centre of the cigar gradually 
moving in if down until it reaches the non-perturbative region at Y ~ bt 2 , it turns out 
that for Y ~ t an abrupt transition occurs and the dominant trajectory switches from 
being centred around if — t (fig. [I]) to being centred around t' = i, where i > is in the 
non-perturbative region (fig. ^b). For such a transition, which we named 'tunneling', 
there is no passage through any intermediate 'banana' configuration (fig. ||]a). 

In this article we wish to explain this phenomenon in more detail than was done 
previously and for the case of the full BFKL equation with running coupling. In sec- 
tion |2| we present analytical calculations which give the estimate of the nonperturbative 
Pomeron, while in section [3] we illustrate the phenomenon with a number of figures 
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resulting from the numerical solution of the BFKL equation. 



2 Analytical considerations 



It is useful to recall the form of the BFKL equation as we shall use it here. We define 
a BFKL Green's function G(Y, t, to) through the following integro-differential equation 



dG(Y,t,t c 
dY 



d 2 Q - n 2, 
— ^a s (lng J 
irq z 



^G{Y, In k'\ t ) - G(Y, t, t Q )Q(k - q) 



k 2 = e l 



k' 



with the initial condition 



G(0,t,t ) = 6(t-t ), 



\k + q\, 
(1) 



(2) 



and with a s = a s Nc/^- All scales are in units of Aqcd- We shall take a s as having its 
one-loop perturbative form, regularised by a cutoff at a non-perturbative scale i: 



a s (t) 



eg - 1) 



12 



(3) 



In some contexts we will use an equation similar to (|1|), but with a s (lng 2 ) replaced by 
a s (t). 

It will be convenient to define some high-energy exponents: cu s , the perturba- 
tive saddle-point exponent, determines the high-energy behaviour of the Green's func- 
tion when it is dominated by scales of order t. It is given by u s = (41n2)a s (t) + 
O(ol). We shall also need the non-perturbative pomeron exponent, tup, which gov- 
erns the asymptotic high-energy behaviour of the Green's function, G(Y) ~ exp(o;py), 
once the problem has become dominated by non-perturbative scales. It is given by 
lujp = (41n2)a s (f) + O (a s (i) 5 ^ 3 ), though it is to be kept in mind that the expansion in 
powers of a^ 3 converges very slowly. 

For the purpose of our arguments we shall take the formal perturbative limit of 
t ~^>t and correspondingly u s u-p. 



2.1 Informal arguments for tunneling 

The mechanism of tunneling is illustrated in figure 0b. To establish when tunneling 
takes place we have to equate the contribution from the tunneling configuration with 
that from the traditional perturbative configuration. Neglecting all prefactors, the 
tunneling configuration gives a contribution 

G tunncl ~e-^/ 2 e-^e-^/ 2 , (4) 

corresponding (reading from left to right) to the 'penalty' for branching from t to i, 
the BFKL evolution at i and the penalty for branching back to t. In contrast, the 
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perturbative contribution is proportional to 

^ ~ e^ Y . (5) 

Equating these two expressions gives 

^tunnel(^) = 7T\ • (6) 

— UJ s [t) 

Perhaps the most important characteristic of this equation is the fact that F t unnei is 
proportional to t rather than t 3 as for the diffusion mechanism. One possible exception 
is the case in which the Pomeron is weak, that is u-p <C 1. Then i should be quite large, 
and including in Eq.(H) the value co s (t) = Xm/ (bt) and the rough estimate wp = XmJ (bt) 
yields yt unne i = (bt/x m )i, so that the perturbative regime is kept up to u s (t)Y values of 
order i. 



2.2 Formal arguments 



For an isolated Pomeron "bound state" , the tunneling probability is estimated on the 
basis of the approximate decomposition 



G(Y;t,t ) = G pe ^Y;t,t )+e^ Y ^p^ 



(7) 



where T Uw (t) is the regular solution for t — > oo of the BFKL equation |24[], satisfying 
the strong-coupling boundary condition also. The decomposition (|7|) separates the 
continuum from the lowest bound state contributions to the gluon's Green functions 
[ 22 1 or, equivalently, the background integral from the leading Regge pole |2l |. The 
normalisation of the bound state factor in front of e UwY is simply obtained by imposing 
that it be a projector by t Q integration, or equivalently from the requirement that 



G(Y;t,t ) = J dt'G(Y;t,t')G(Y-y;t',t ). 



(8) 



Eq. (0) applies as it stands to a symmetrical BFKL kernel, as that with the scale choice 
a s (lng 2 ) in Eq.([I]). If instead the scale a s (t) is chosen, the kernel is asymmetrical, the 
conjugated wave function is F^it) = tJ 7 * (t), and the bound state projectors differ by 
a to factor and by the normalisation integral 



G(Y;t,t ) = G pert (Y;Mo) + e^ Y ^]^l to . 



(9) 



In either case, the large t behaviour of T^^t) is estimated in a semiclassical approach 
in which ba s (t), 1/Y and also bu-p are taken to be small parameters (weak Pomeron, or 
b - expansion p2[ ). In such a case the 7 - representation of is valid, and takes the 
form 



1/2+ioo 7 

' 07 (7-1/2)*- 



1/2-ioo 



2m 



bujp 



-27r X '[p(b^t)} 



exp 



bid 



II' 



bu)-pt 



p(x) dx 



(10) 
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where the saddle point condition at 7 = | — p 



X\ 2 



P 



ax( 7 ) 



buj-pt 



(11) 



defines the semiclassical "momentum" p(bu-pt). The effective eigenvalue function in 
Eq . (1TT| ) is just the leading characteristic function Xo(l) f° r the a s (t) choice, and differs 
by some corrections involving the parameter wp in the symmetric scale choice a s (lng 2 ) 
which essentially lead to a rescaling of the (right) regular solution by a constant factor 
~ VVXm (see Appendix) . 

The Pomeron suppression arising from Eq.flTOD has a universal form ~ w^w^M*))) ^ 
where the function g(a s ) (g(0) 7^ 0) is easily found from Eq. flTTD , [2^|. On the other 
hand the perturbative part, corresponding to the continuum, is dominated by the saddle 
point exponent uJ s (t) 



Gpertit, tf)] Y) 



1 



exp ( uj s \ j Y + corrections^ , (12) 



where the diffusion corrections are due to the running coupling ]T8|, [19], ^ • 

From Eqs.(|7|j9|fr0|) we estimate the t-dependence of the Pomeron contributions for 
the choices d s (k 2 ) (a s (q 2 )) of the running coupling scale respectively 



G 



(k) 
P 



exp 



G 



(<i) 
p 



2-nbuspt 
27r(Mpt) 2 Vx 7 



-t+— (logMpt + 1 + V(l) + O ( 77 * 3 



p 



exp 



-t + — logMpt + 1 + ^(1) + O 



(bupty 



Notice that they differ both in the t-dependence, because of the t factor on Eq.( 
for the t-dependent normalisation integrals 



(13) 
and 



iV«(t)=(J_ dtt^Jt) 



(14) 



The latter quantities cannot be estimated in the semiclassical approximation, which 
breaks down at the turning points, and they need to be evaluated numerically. We can 
roughly estimate their size on the basis of the Airy diffusion model []24j obtained by the 
quadratic expansion of the characteristic function 



XM = Xm + 2Xm(7-7m) 2 + 



(15) 
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This approximation is expected to be a reliable one for the Pomeron dynamics provided 
the fluctuations ((7 — 7 m ) 2 ) ~ "^r- are small compared to 4!%^/ (2xm ), and t varies in 
a range for which 



2bu 



y" 

Am 



Xr, 



buj-p 



(16) 



is a parameter of order unity. In this case the BFKL equation with the choice a s (k 2 



reduces, by Eq.([L5]) to a second-order differential equation in t [24|, and the cut-off 
condition for a s (t) at t — i @ implies that the regular solution for t — > — 00 must 
vanish at i , and for t around i takes the following form 



y" N 5 

Am 



2ku 



7T 



Bi[^(*)] - ?f#Ai[£^)] 



Ai[&(*)]" 



while the regular solution for t — > 00 becomes 



v 

At, 



2ku 



^(t) = Ai [£,(*)] 



with relative normalisation fixed by the Wronskian 



Am 



(17) 



18) 



(19) 



Using now the Airy form for the regular solution fllB]) we can evaluate the normali- 
sation integrals (|T^) to get 



N ik \i) 



Xm ( 2buj w ,\ 
buJTP V Xr 



Xi 



eo[Ai'(-^ )] 



) 



2 3 



5() 



de[Ai(0] 2 + / dee[Ai(0] 



6) 



^ s (t) + 



o; s (f) - 



d f Al (0 =^0^tx 



[Ai (~£o)] 2 , (20) 



where we have used the Airy identities [ 25| and the assumption that the Pomeron 
intercept is given by the rightmost zero of the Airy Ai(£) function 



(2buj 



TP 



\ y" 

\ Ar, 



Xr. 



bujp 



-fo ^ -2.3381 



Ai(-£o 







(21) 



The additional factor \m which multiplies N^(t) comes from the rescaling of the 
(right) regular eigenfunctions due to a scale change ^(t)y^ — J~$£(t). 

The expressions ( |T3"D supplemented by the normalisation factors ( |211D can be com- 
pared with the naive estimate (H) by rewriting them in the form 



G 



(k) 
ip 



G 



(9) 
p 



a s (i) exp 
«s(t) exp 



ocjp r 



(22) 



where one sees the explicit a a (t), (a s (t) 2 ) dependence for a s {k 2 ) and a s {q 2 ) scale choices 
respectively. The log ~ in the exponent represents a collinear enhancement due to the 
evolution from scale t to scale i, and the i - dependent factors 



exp(A< fc >) 
exp(A (9) ) 



1 



27r£ [Ai'(-£ )] 2 + fw s (t) cup 



exp 



6u; 



;i+V(i)+iog6wpt)-t 



IP 



[^ s (t) - CUp] 



27r£o[Ai'(-£ )] 2 w 8 (t)u; 



exp 



ii' 



;i + ^(i)+iog&wjp£) 



(23) 



come from the normalisation factors (|20|) and of the semiclassical prefactors. By then 
comparing the result (22) with the perturbative Green's function in Eq.(|12"D we obtain 
the estimate of the tunneling transition rapidity in both scale choice cases 



(u F -uj s (t))Y t 



(k) 

tunnel 



(utp - uj a (t)) Y t 



(?) 

tunnel 



2 t 1 
t-t- - — log = - log a s (t) - - log(27rx^Q! s (t)y tunne i) - 
ocjp t 2 

2 t 1 
t-t- log = - 2 logout) - - log(27Tx^a s (i)ytunnci) 



A«(f) , 

-AW(t), 
(24) 



2.3 Limits of the tunneling picture 

We have just seen that the Pomeron regime sets in, according to Eqs.(|6]) and ( pi"D at 
u s (t)Y values of order t. This is to be compared with the validity boundary of the hard 
Pomeron behaviour |22] 1 < uo s (t)Y < t/y/D, where D = ^j— is the diffusion coefficient. 
Thus, for phenomenological values of t and t there appears to be little room for the 
hard Pomeron (and its diffusion corrections) to be observable, at least for models (such 
as the leading BFKL equation) in which the hard Pomeron and the nonperturbative 
Pomeron are quite strong. The Pomeron contamination in such cases, is a serious 
problem even for a theoretical determination of perturbative diffusion corrections to 
the hard Pomeron. In fact, one has to use the ^-expansion |22] -6^0 limit with 
a s (t) and a s (t) fixed — in order to identify such corrections at the numerical level. 

The situation may change, if subleading and unitarity corrections are taken into 
account. Subleading contributions to the kernel are large |J and need to be resummed 
0, ||, U : they provide eventually a sizeable decrease of the hard Pomeron intercept and 
even more of the diffusion coefficient. This effect slows down both the nonperturbative 
effects and diffusion corrections, thus increasing the range of validity of the perturba- 
tive predictions. Furthermore, unitarity effects are expected to be important for large 
densities and more so in the strong coupling region. All this goes in the direction of a 
weaker asymptotic Pomeron — as is the actual physical Pomeron of soft physics [f26 



thus increasing the effective t value, and the window of validity of the hard Pomeron 
regime.^] 

2 Indeed if the soft Pomeron has weaker energy dependence than the hard pomeron, then tunneling 
stops being a favoured contribution altogether since the non-perturbative Pomeron cannot overtake 
the hard pomeron before the oscillatory regime is reached. 
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It may thus be that in realistic small- a; models |14| with a weaker Pomeron, diffusion 
corrections to the hard Pomeron may set in before tunneling, so that the transition to 
the Pomeron regime is smoother. In any case, before diffusion corrections lead to a 
decreasing cross-section and to the asymptotic oscillatory regime [21] the Pomeron is 
expected to take over. 



3 Pictorial results 
3.1 Cigars 

Perhaps the most dramatic illustration of the phenomenon of tunneling can be obtained 
through pictorial representations analogous to those of figures [I] and 0, but obtained by 
exact numerical solution of the BFKL equation. 

This involves studying the running-coupling dynamics of BFKL evolution and in 
particular the sensitivity of the Green's function G(Y, t, t) to evolution at some inter- 
mediate point y,t'. This can be done by examining: 

/ (y . g .«.0 = g^g^ , (25) 
which, by virtue of the linearity of the BFKL equation, is normalised as follows 

dt'f(Y,y,t,t') = l. (26) 



From the point of view of graphical representation this quantity suffers from the problem 
that for y ~ or y ~ Y and t' ~ t one is sensitive to the 5-function initial condition for 
G, eq. (@). To smooth-out this 5-f unction we examine instead a smeared version of /: 

¥(Y _ fdtodhe-^/^-tf/iGjy^^t'MY- y,t',t) 

I[ ' V " ' ] Jdt dt x e -^-t) 2 /2-(t,-t)V2 G(Y, t , h) 1 } 

Contour plots of f(Y, y, t, t') as a function of y and t' are shown in figure ||] for various 
values of Y. We have chosen a value for t that is physically unrealistic, t = 30. The 
reason for this choice is that it enhances the separation between the perturbative and 
non-perturbative regions, allowing us to illustrate more clearly the transition that takes 
place.0 The coupling is regularised in the infrared, eq. (||), by a cutoff at i = 4. 

At the lowest value of Y in figure 13 the shape of the cigar is largely determined 



by the smearing of the initial condition, eq. (27). At Y = 50 one can see a significant 
amount of diffusion at central values of y, as well as a slight asymmetry favouring lower 
t' values, caused by the running of the coupling. 

At Y = 60 the tunneling transition starts to take place: at low values of t' ~ t, 
an enhanced region appears which is disconnected from the main region of the cigar 

3 For more phenomenologically relevant plots we note that in addition to taking more reasonable 
values for t it would also be necessary to include the NLL corrections and collinear resummation. 
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Figure 3: Contour plots for f(Y,y,t,t'), illustrating different stages of the evolution: 
Y = 30 and 50 illustrates standard 'cigar' type plots, Y = 60-70 show the point 
where tunneling begins to play a role (BFKL 'aliens'), while for Y > 80 tunneling has 
become the dominant contribution (BFKL 'bowls'). MPEG animations of the transition 
are available from http : / /www . lpthe . jussieu . f r/~salam/tunneling/ which includes 
additionally results for other combinations of t and t. 
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- it has not appeared through diffusion, but through tunneling. As Y increases the 
perturbative region of if ~ t becomes progressively less important, eventually being 
completely supplanted by the non-perturbative region t' ~ i. In this transition stage 
we take the liberty of referring to the contour plots as 'BFKL aliens'. 

Beyond Y = 75, the evolution takes place entirely in the non-perturbative region, 
except in the wings, y ~ and y ~ Y, which connect the initial conditions at t to the 
infrared region. We refer to these contour plots as 'bowls'. 

We note that the position in Y of the tunneling transition is quite consistent with 
the arguments of section |^. For our particular non-perturbative conditions, the non- 
perturbative power is cup ~ 0.32; substituting this into eq. (p4[) gives Y tunnc \ ~ 68. 



3.2 Ducks 

A convenient way of visualising the tunneling transition within a single plot is to con- 
sider f(Y, Y/2, t, t') as a function of Y and if , i.e. a contour plot built up of cross sections 
through the centres of the plots of figure ||], for a range of Y values. This is shown in 
the upper plot of figure At low Y we see just the width of the initial condition, which 
gradually diffuses out while remaining centred at if ~ t. At Y ~ 60 a second region 
appears 'out of nowhere' at if c± i and the original cigar region disappears completely 
just beyond. This, once again, is the tunneling transition. We refer to the resulting 
figure as a 'BFKL duck'. 

In the lower plot of figure [| we show two quantities as a function of Y (on the same 
Y scale as the upper plot): the effective exponent u e fi, defined as 

^ = 4f ln l dt ° dk e~ {t °- t)2/2 - {h - t)2/2 G(Y, t , *0 , (28) 

and the position of the maximum in t' of f(Y,Y/2,t,t'), which we refer to as t max - At 
the point where the tunneling transition occurs there is a jump in the value of t max 
as it switches from the perturbative to the non-perturbative region. Simultaneously 
there is a rapid change in u e g as it goes from the perturbative exponent ~ uj s to the 
non-perturbative one u-p. 

We can use the point where the position of the maximum goes from being above to 
below if — |(t + 1) as the definition of where tunneling occurs, Itunnei (in what follows, 
since we will no longer be faced with the problem of visualising an initial 5-function, we 
will actually consider f(Y,Y/2,t,t') rather than /, though the difference in practice is 
minimal). Alternatively one can define the tunneling transition as occurring when the 
integrals of f(Y, Y/2, t, t') above and below |(t + t) are equal, 

/(t+t)/2 ^oo 
dt'f(Y, Y/2, t, = / dt'f(Y, Y/2, t, t') . (29) 
-oo J(t+t)/2 

These two definitions for Ytunnei have been used to deduce, from our numerical solutions, 
the points with errors bars in fig. [5] (the lower error bar corresponds to the first definition, 
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Figure 4: The upper figure is a contour plot of f(Y, Y/2, t, t'), corresponding to a cross 
section at Y/2 of 'cigar'-type plots in fig. |3|, as a function of Y. We refer to the resulting 
figure as a 'duck'. The lower plot shows the effective exponent, u cS , and the position 
of the maximum of /, t max) as a function of Y. 
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while the upper error bar corresponds to the second one), where we plot Itunnei as a 
function of t, for two values of i. They are to be compared to the smooth lines, which 
correspond to the analytical predictions both in a s (t) and a s (lng 2 ) cases. We have 
tried the direct evaluation of the Gp from Eqs. ([?[), ([J) by numerically calculating the 
7-representation integrals (solid lines) and the analytical approximation based on the 
Airy model Eq. (|24D (dashed lines). 
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Figure 5: The position in Y of the tunneling transition, Ytunnei as a function of t for 
two different non-perturbative cutoffs i (points). The points represent the numerical 
'measurements' of yt unnel based on the criteria discussed in the text. They are compared 
to the prediction of Eq. (p4"|), dashed lines and the calculation which involves the direct 
numerical evaluation of the 7-representation in the normalisation integrals Eq. (|l4l), 
solid lines. Left-hand plot: the calculation for a s (lng 2 ) scale choice, right-hand plot: 
a s (ln k 2 ). 



One sees a remarkable agreement between the 'numerically measured' and analytical 
results, which would seem to be strong confirmation of even the fine details of our 
picture. We note in particular the roughly linear dependence of Y tnnnc \ on t and the 
fact that its slope increases for larger i as a consequence of the smaller ujp and hence 
smaller denominator in eq. (|6]). 



4 Conclusions 

In this paper we have presented a detailed investigation of the abrupt tunneling tran- 
sition from the perturbative regime to the non-perturbative 'Pomeron' regime for the 
case of leading-log BFKL evolution with a running coupling. 

The tunneling transition takes place when u s Y ~ i, with exp(f) the effective non- 
perturbative scale of the problem (in units of A.q CD ) . This is parametrically much earlier 
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than estimates of the maximum perturbatively accessible Y value that are based on 
diffusion arguments. Furthermore, in contrast to expectations from diffusion arguments, 
the transition does not involve evolution at scales intermediate between those of the 
hard probe and the non-perturbative region, as shown in figs. |3] and 

Since the purpose of this paper has been to give an illustration of the effect of 
tunneling, we have restricted ourselves to leading-logarithmic evolution and considered 
only very extreme scales for the hard-probe, where the effect is more dramatic because 
of the large ratio up/w a (t), or equivalently t/t. 

In a more phenomenologically relevant context there remain a number of other 
important issues: the separation between tunneling and diffusion may be less distinct, 
both because the accessible range of t is smaller, and because for a given ratio of t/i 
higher-order corrections tend to reduce ujp/u s (t). Furthermore unitarity and saturation 
corrections may actually cause the effective value of uo-p to be smaller than io s which 
could eliminate tunneling altogether. Developing a better understanding of these and 
related questions remains one of the main challenges of small- x physics. 
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Appendix A: Pomeron weight with cx s (\nq 2 ) 



Taking a s (\nq 2 ) instead of a s (t), as in Eq.([l|) means adding a term — §(xo + Xo) t° the 
NL kernel which emerges due to the evolution from scale t to scale q 2 . Correspondingly, 
in the 00 - expansion formulation || ^§ the effective characteristic function becomes 
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As a consequence, the large t-anomalous dimension stays the same 
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It implies that approximately the eigenfunctions JuQ(t) get rescaled by the additional 
constant normalisation factor so that ^m^i* — Fuw- 
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